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I. 


INTRODUCTION 


Phase-locked  loops  are  often  employed  to  derive  synchronous  local  carriers 
for  coherent  communication  systems.  Because  the  stabilities  of  the  various  system 
oscillators  can  limit  performance,  the  design  engineer  is  often  called  upon  to  specify 
oscillators  so  that  they  will  cause  less  than  a  given  amount  of  performance  degradation. 

In  this  paper  we  shall  examine  some  effects  of  short-term  oscillator  instabilities  upon 
the  second  order  phase-locked  loop.  We  also  show  how  the  magnitude  of  these  effects 
may  be  predicted  from  simple  measurements  or  from  conventional  RMS  frequency 
deviation  specifications.  Very  little,  if  any,  of  the  content  is  new  as  such;  the  utility 
of  this  paper  resides  in  demonstrating  the  interrelations  among  the  results  of  various 
workers  and  providing  guidelines  for  the  short-term  stability  specifications  of 
oscillators  for  coherent  communications  systems. 

II.  COMMUNICATION  SYSTEM  MODEL 

To  provide  a  vehicle  for  discussion,  we  assume  that  the  system  under  study 
may  be  reduced  to  the  simple  block  diagram  of  Fig.  1.  The  unmodulated  transmitted 
signal  suffers  a  fixed  attenuation,  is  delayed  by  a  (possibly  time -varying)  phase  shift 
and  is  received  in  the  presence  of  additive  Gaussian  noise  n(t).  The  spectrum  of 
the  noise  is  assumed  flat,  with  real  density  Nq  watts/Hz  across  the  bandwidth  of  the  RF 
filter.  No  fading  is  allowed  —  the  desired  received  signal  has  constant  power  P^. 

The  word  oscillator  is  used  to  denote  any  set  of  frequency  synthesis  equipment. 
Thus  the  box  labeled  transmitter  oscillator  may  correspond  to  an  actual  interconnection 
of  several  oscillators,  amplifiers,  mixers,  filters,  etc.  ,  whose  nominal  output  frequency 
0)^  is  the  desired  carrier  frequency.  Similarly,  the  receiver  equivalent  voltage  con¬ 
trolled  oscillator  (VCO)  would  ordinarily  be  realized  by  one  or  more  oscillators,  mixers, 
filters,  etc. 

The  received  signal  is  written  as  /2P~  cos  [<j  t  +0^  (t)  +^>(t)]  +  n  (t). 

0^  (t)  represents  phase  variations  of  the  transmitted  signal  while  cp( t)  represents 
channel -induced  phase  fluctuations. 

III.  EQUIVALENT  OSCILLATOR  MODEL 

The  short-term  phase  fluctuations  observed  in  the  output  of  a  quartz  crystal 

12  3 

frequency  standard  are  attributed  to  three  major  effects.  ’  ’  One  type  of  fluctuations 
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result  from  perturbations  of  the  oscillator  by  thermal  and  shot  noise  generated  in  the 
crystal  and  associated  oscillator  circuitry.  A  flat  power  spectral  density  of  frequency 
fluctuations  is  associated  with  this  noise. 

A  second  type  of  fluctuations  is  attributed  to  crystal  unit  and  circuit  parameter 
changes.  The  resulting  power  spectral  density  of  frequency  fluctuations  appears  to  have 
an  f  1  characteristic.  (Of  course,  the  true  power  spectrum  cannot  behave  as  f  ^ 
all  the  way  down  to  f  =  0,  because  this  would  imply  infinite  power  in  the  process  and, 
therefore,  arbitrary  long-term  frequency  excursions.  In  the  following  analysis  the 
spectrum  is  always  weighted  by  a  function  having  a  zero  at  f  =  0. )  This  type  of  noise 
is  often  referred  to  as  flicker  noise. 

The  third  type  of  fluctuations  is  produced  by  additive  noise.  To  obtain  long-term 

stability,  oscillators  usually  operate  at  low  power  levels  and  are  followed  by  several 

amplifier  stages  to  achieve  a  reasonable  output  power.  These  amplifiers  add  noise 

(assumed  Gaussian)  to  the  signal.  Often  a  narrowband  filter  follows  the  amplifiers  to 

reduce  this  noise.  If  the  noise  power  is  small  compared  to  the  desired  signal  power, 

half  the  noise  appears  as  amplitude  modulation  sidebands  and  the  other  half  as  phase 

2 

modulation  sidebands.  The  resulting  power  spectral  density  of  the  frequency  fluctua¬ 
tions  depends  on  the  particular  clean-up  filter  used. 

Temperature  variations,  automatic  gain  control  perturbations,  power  supply 
and  loading  changes,  and  vibration  can  all  contribute  additional  phase  fluctuations  to 
oscillator  outputs. 

Perfect  frequency  multiplication  by  a  factor  of  N  will  multiply  phase  variations 
by  the  same  factor.  Aside  from  this,  the  primary  degradation  of  the  short-term 
stability  of  the  oscillator  introduced  by  well -designed  frequency  synthesis  equipment 
should  be  additive  noise.  Therefore,  the  general  form  of  the  spectrum  of  the  phase  of 
the  frequency  synthesis  equipment's  output  may  be  taken  to  be  the  same  as  that  of  a 
crystal  oscillator. 

Now  if  we  let  the  oscillator  (in  the  general  sense)  output  be  cos  [  cj^  t  +  0  (t)]  , 
and  if  we  assume  an  output  filter  whose  equivalent  positive  frequency  transfer  function 
is  given  by 

j(w  -  £0c)  +  Uj 
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then  the  assumed  power  spectral  density  for  the  frequency  fluctuations  is  given  by 


sew'Ni+]q+N3 
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N.  is  the  equivalent  white  noise  density  of  the  perturbing  noise,  N9  is  the 

1  -1  1 
constant  describing  the  intensity  of  the  f  noise,  and  is  the  effective  unfiltered 

additive  noise  density.  The  equivalent  phase  noise  density  is  given  by 
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The  assumption  used  above  that  the  additive  noise  power  is  much  less  than  the 
signal  power  corresponds  to  satisfaction  of  the  inequality  ^  «  1. 

IV.  PHASE -LOCKED  LOOP  MODEL 

The  following  discussion  assumes  some  familiarity  with  phase-locked  loops. 

More  detail  may  be  obtained  from  Ref.  4. 

In  normal  operation,  the  phase -locked  loop  error  0  will  be  much  less  than  1 

radian,  so  that  a  linearized  analysis  is  sufficiently  accurate.  For  random  disturbances, 

2 

the  variance  of  the  phase  error  a  9  provides  a  convenient  measure  of  performance; 

€  2 

operational  thresholds  are  often  specified  by  maximum  allowed  values  of  Oq  .  Clearly 

2  £ 
the  part  of  Qq  ~  due  to  oscillator  instabilities  should  be  much  less  than  that  due  to  the 

additive  noise  and  channel  fluctuations  at  the  minimum  received  signal -to -noise  ratio, 

if  the  oscillators  are  not  to  degrade  system  performance. 

The  loop  filter  transfer  function  is  assumed  to  be  of  the  form  K(1  -I-  a/s).  The 

loop  damping  ratio  £  and  natural  frequency  to  are  related  to  these  parameters  by 
2  n 

2£  =  A  K  and  =  a  A  K.  The  effective  gain  A  is  given  by  A  =v/^P^  • 

The  equivalent  linearized  model  for  the  system  in  Fig.  1  is  given  in  Fig.  2. 

(The  details  of  the  derivation  of  the  equivalent  phase-locked  loop  model  are  given  in 
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Ref.  4.)  The  effective  noise  n'(t)  is  modeled  as  white  Gaussian  noise  with  two-sided 

N0 

spectral  density  -j- .  The  subscripts  T  and  R  denote  transmitter  and  receiver, 
respectively. 

The  oscillator  perturbing  noises  are  represented  as  n^.  (t)  and  n^R  (t).  The 

noises  n^  (t)  and  n  (t)  represent  oscillator  flicker  noise,  while  n  (t)  and 

n^R  (t)  represent  the  effects  of  additive  oscillator  noise.  The  noises  n^  (t)  and 

n.„  (t)  are  included  to  describe  miscellaneous  effects  such  as  vibration.  Note  that  all 
4R  ' 

of  these  terms  represent  the  frequency  fluctuations  due  to  the  corresponding  physical 
processes. 

The  spectra  associated  with  the  first  three  oscillator  noises  are  given  by 
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Writing  out  the  operator  equations  describing  Fig.  2  we  find 
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The  first  term  describes  the  effects  of  the  additive  noise,  while  the  second  shows 
the  effects  of  channel  phase  fluctuations.  The  third  term  shows  the  effect  of  the  oscillator 
noises.  Note  that  the  transmitter  and  receiver  oscillator  noises,  and  channel  phase 
perturbations,  affect  0  in  the  same  manner. 

V.  LOOP  PHASE  ERROR  VARIANCE 

Since  the  various  noises  are  assumed  independent,  the  loop  error  variance  may 
be  written  as 


T  2  ;  2 

\  ih  ai 


where  a,  is  due  to  received  additive  noise,  cr9  is  due  to  random  channel  phase 
1  2  2  L 

fluctuations,  and  a g  and  a  £  are  the  contributions  of  the  transmitter  and  receiver 
oscillators,  respectively.  By  direct  calculation, 
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where  S^  (cj) 


is  the  power  spectral  density  of  the  channel  phase  variations. 
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where 
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S  (u)  is  the  power  spectral  density  of  the  noise  due  to  miscellaneous  effects. 
n4T 

The  function  f(£)  is  given  by 


*<  1 


f(0=  <  1 


i  =  i 

s>i 


The  expression  for  is  the  same  as  that  for  cr^  with  R  substituted  for  T. 
A  common  choice  for  the  loop  damping  ratio  is  £  =  ,  for  which  we  get 


6 


We  observe  that  for  sufficiently  small  loop  natural  frequency  u  and  miscel- 

2  n  2 

laneous  transmitter  oscillator  noise  variance  O'  the  term  .25  N-_A>  (due  to 

2  vT  2T  n 

flicker  noise)  dominates  while  for  sufficiently  large  u  either  the  term  containing 

o  n 

(due  to  perturbation  noise)  or  the  term  containing  (due  to  additive  noise) 

will  dominate.  In  many  applications  the  term  corresponding  to  the  perturbation  noise 
12  5 

can  be  neglected.  ’  ’ 

For  typical  carrier  tracking  loops,  the  ratio  of  the  loop  natural  frequency  to  the 

post-oscillator  noise  filter  (u  /u^kis  very  small,  so  that  the  third  term  in  the 
2  n  IT  N3t 

expression  for  cr^  becomes  — 1 — ^ - ,  which  is  the  total  phase  noise  power 

due  to  the  additive  Gaussian  noise. 

2 

The  above  discussion  also  applies  to  a ^  due  to  the  receiver  oscillator. 

However,  in  drawing  the  block  diagram  of  Fig.  2  we  have  assumed  that  the  phase  shift 
characteristic  of  the  filter  following  the  local  oscillator  could  be  ignored.  This  implies 
(u^/u  ^)  «  1.  If  this  condition  is  not  satisfied,  the  effect  of  the  filter  upon  loop  stability 
and  tracking  dynamics  must  be  examined. 


VI.  MEASUREMENT  OF  cr  *  AND  a £ 

2  2 

The  maximum  allowed  values  of  cr^  and  cr^  ,  together  with  the  loop  damping 

ratio  £  anc*  natural  frequency  cj  ,  define  the  necessary  transmitter  and  receiver 

n  2  2 
short-term  stability  specifications.  A  typical  requirement  might  be  a 3  <  .005  (rad)  , 

under  given  environmental  conditions,  which  insures  that  the  loop  threshold  performance 

is  controlled  by  the  received  signal -to -noise  ratio  and  received  signal  dynamics.  This 

2  2 

section  of  the  paper  describes  several  techniques  for  the  measurement  of  cr^  and  O ^  . 

2  2 

A.  Direct  Measurement  of  a ^ 


If  the  transmitter  and  receiver  are  available  in  the  same  place,  and  if  the 

2  2 

receiver  phase-locked  loop  gain  may  be  temporarily  modified,  a ^  maY 

measured  as  shown  in  Fig.  3.  The  signal  e(t)  is  the  phase-locked  loop  phase  detector 
output.  The  receiver  phase -locked  loop  gain  must  ordinarily  be  changed  to  correct  for 
the  variation  of  £  and  with  the  received  signal  strength. 


7 


Of  the  various  types  of  averaging  techniques  which  might  be  used  in  the  con¬ 
figuration  of  Fig.  3,  a  particularly  convenient  one  employs  a  low  pass  RC  filter  with 

transfer  function  H(s)  =  1  *  .  Then  the  expected  value  4  of  the  averager  output 

2  -L  1  1  s  2  2 

y(t)  is  E[e^  (t)],  which  is  Og  +  CJ4  •  Of  course  since  y(t)  is  a  sample  function  of 
a  random  process,  it  can  (and  generally  will)  vary  with  time. 

The  time  constant  T  may  be  selected  so  that  |y(t)  —  n\  is  less  than  some  desired 
value  for  a  specified  percentage  of  the  time.  The  primary  difficulty  in  giving  an  a 
priori  specification  of  T  is,  of  course,  that  the  spectrum  of  e(t)  must  be  known.  In 
practice  a  T  on  the  order  of  several  ^  seconds  should  be  sufficient.  A  detailed 
discussion  of  bounds  on  the  variance  of  y(t)  is  given  in  Ref.  6. 

2  2 

B.  Direct  Measurement  of  or  q 4 

If  the  transmitter  or  receiver  oscillator  is  replaced  by  a  frequency  source  known 

2  2 

to  be  very  good  the  arrangement  of  Fig.  3  may  be  used  to  obtain  Qo  (or  aA  )  directly. 

J  2 * 2 *  4 * * 7  2 

If  the  replacement  source  is  not  particularly  good,  the  true  value  of  0^  or  04  may  be 

upper  bounded  by  assigning  all  of  the  measured  phase  jitter  power  to  the  system 

2  2 

oscillator.  If  the  oscillators  are  similar,  one-half  of  the  total  variance  +  04  is 
sometimes  attributed  to  each  oscillator. 

Often  one  must  determine  if  the  stability  of  a  given  oscillator  is  adequate  without 

direct  access  to  either  the  complete  receiver  or  transmitter.  Then  one  approach  would 

be  to  construct  a  phase-locked  loop  whose  error  behavior  is  designed  to  be  identical  to 

7 

that  of  the  receiver,  as  indicated  in  Fig.  4. 

However,  because  the  transfer  functions  from  the  transmitter  and  receiver 

oscillator  noise  sources  are  the  same,  the  oscillators  under  test  may  be  interchanged 

with  the  reference  oscillator  of  Fig.  4  to  give  the  alternate  test  configurations  shown 

in  Fig.  5.  From  the  second  configuration  of  Fig.  4  and  the  first  configuration  of 

2  2 

Fig.  5  we  note  that  only  one  test  setup  is  needed  to  measure  both  a 3  and  04  . 

C.  Discriminator  Test 

Often  for  one  reason  or  another,  a  test  configuration  employing  a  phase -locked 
loop  is  inconvenient.  If  a  discriminator  is  available,  the  alternate  technique  illustrated 
in  Fig.  6  may  be  used.  The  filter  transfer  function 
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where  G  is  a  scaling  factor,  is  selected  so  that  the  linearized  transfer  function  from 
the  equivalent  oscillator  noises  to  the  equivalent  error  voltage  e(t)  is  identical  to  that 
for  the  original  receiver. 

VII.  RMS  FREQUENCY  DEVIATION  SPECIFICATION 

Short-term  oscillator  stability  is  often  specified  by  a  measured  RMS  frequency 

deviation  ^  for  a  given  averaging  time  r.  In  this  section,  we  wish  to  point  out  how 

the  parameters  N.  ,  N9,  and  N~  characterizing  the  short-term  stability  of  a  particular 
1  L  6  Af 

oscillator  may  be  determined  from  given  —  versus  r  data. 

Af  1 

The  quantity  —  is  related  to  oscillator  behavior  in  the  following  manner.  Let 
the  oscillator  output  signal  be  cos  [cjq  t  +0  (t)]  where  0  (t)  represents  phase  fluc¬ 
tuations.  The  average  oscillator  frequency  departure  <0>  of  a  time  interval  r 
seconds  long  centered  on  time  t  is  defined  by 

t  +  r/2 

<0>=7  "  0(g)  d /u=^ 

t,T  T  t  —  t/  T 


where  A  0  is  the  difference  0  (t  +  r/2)  —  0(t  —  r/2).  Let  0(j  r  +  t)  =  0^  and  A0^  = 

0  —  0^  j  .  The  measured  variance  of  N  adjacent  samples  of  the  average  frequency 
departure  is  determined  from 
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is  included  to  eliminate  the  effect  of  a  constant  drift.  The  average  value  of  <j~  (N,  t) 
is  related  to  the  power  spectral  density  of  the  phase  fluctuations  by 
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Then  the  RMS  frequency  deviation  is  given  by 


Af  =  Jo2  (N,  r) 

f  0 

o 


Note  that  several  experiments  should  be  performed  to  get  a  (N,  t)  .  See  Ref.  2  and 
Ref.  8  for  more  details. 

Now  if  we  assume  the  form  of  the  power  spectral  density  of  the  frequency 
fluctuations  is 


2  2 

N„  No w  X 
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as  described  above,  we  find 
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(N,  r) 


-  +  K  N0  +  - rj- 

T  2  2t 


T 


where  K  =  .66  +  .318  In  (N/2)  ;K  =  1.90  for  N  =  100. 

A  typical  asymptotic  plot  of  y-  versus  r  is  given  in  Fig.  7.  The  dominant 
contributions  are  labeled  for  various  ranges  of  r.  For  many  frequency  standards  it 

rr  2 

appears  that  the  effect  of  the  perturbing  noise  on  —  is  negligible.  Also,  measure- 
ments  on  VCO  noise  characteristics  ’  indicate  tha?  the  perturbation  noise  gives 
negligible  contributions  to  narrow-band  phase-locked  loop  error  variances  compared 
with  the  contributions  due  to  flicker  noise. 

The  expression  for  a 2  may  be  applied  to  several  measurements  of  y  for 
various  r  and  N  to  provide  estimates  of  the  characteristic  parameters  N^,  N2.  and 

No*  The  contribution  of  each  term  to  0.  for  given  loop  parameters  may  then  be 

J  4  2  2  Af 

determined.  Conversely,  for  a  given  maximum  cr^  or  a^,  bounds  on  y  as  a 

function  of  t  may  be  deduced  from  the  above  equation. 

VIII.  AN  EXAMPLE 

A  phase-locked  carrier  tracking  loop  is  to  be  employed  in  a  telemetry  receiver 
operating  at  250  MHz.  All  local  oscillator  signals  are  derived  by  frequency  multipli¬ 
cation  from  a  single  5  MHz  VCO.  From  dynamic  tracking  and  signal -to -noise  ratio 
considerations  the  desired  loop  parameters  are  £=  .7  and  =  20. 

Because  of  the  way  in  which  the  oscillator  signals  are  related,  the  equivalent 
local  oscillator  signal  is  just  the  VCO  output  multiplied  in  frequency  by  50. 

The  proposed  5  MHz  VCO  has  the  following  performance  for  N  =  100: 

t,  seconds 
1 

.  1 


Af 

f 

10 


-9 


2  x  10 


-9 


11 


.01 


10 


-8 


10 


-7 


.001 


The  1  second  measurement  should  be  dominated  by  flicker  noise,  while  the 
1  msec  measurement  is  dominated  by  additive  noise.  Therefore,  at  5  MHz 


Af  _  a_  ^ 
f  u 

o 


Af 

f 


for  t  =  1  sec 


for  t  =  1  msec 


(The  post -oscillator  filter  has  a  bandwidth  much  greater  than  1  kHz. ) 
Solving:  1.90  N2  =  77^  x  10  ^  and 


N3U1  2  -6 

- * -  =77  X  10 


The  equivalent  parameters  at  250  MHz  are  given  by 


1.  90  N2  =  (50) 2  x  77 2  x  10 '4  =  Q£) 


N3U1  Sir',2  ,n-2 

2 - KYy  X  10 


Since  perturbation  noise  is  assumed  negligible,  we  have 


Ni 

«  1.90  N2 
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which  implies 


Ni  «  Gy 


The  expression  for  becomes 


2  -2'  .4  ' 

CT4  -  1.75x10  +  6.  25  x  1U  4  N2  H — r> — - 


2.  46  x  10 


-2 


N  Q 

The  contribution  — ^ —  due  to  additive  noise  is  much  greater  than  that  from 

either  of  the  other  terms.  Assuming  all  the  other  system  oscillators  are  much  better, 
2 

the  equation  for  (jQ  shows  that  if 


c 


if)  -106 

0  o 


2  2 

is  the  received  signal -to -noise  ratio  required  to  give  Og  =•  1  (rad)  with  ideal  oscillators, 

then  about  1.  33  x  (P  /N  )  =  141  (or  an  additional  1.  2l  db)  is  actually  required  to 
2  o 

obtain  Oq  =  .  1  with  the  given  oscillator  and  specified  loop  parameters. 
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Fig.  1.  Simplified  communication  system. 
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PHASE  -  LOCKED  RECEIVER 

Fig.  2.  Linearized  model  of  simplified  communication  system. 
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Fig.  3. 


Direct  measurement  of  a 
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Fig.  5. 


Alternate  measurements  of 
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Fig.  6. 


Discriminator  measurement  of 
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Fig.  7.  Dominant  effects  upon  short  term  stability.  (  From  Ref.  1) 
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